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1. Introduction 

Fuchsian differential equations (FDE), i. e., homogeneous differential 
equations with single-valued analytic coefficient functions, play a rel- 
evant role in many areas of Physics, Engineering and Statistics. A 
large class of the special functions appearing in these areas are but 
solutions of second order FDE. The most studied of these equations 
is the Riemann hypergeometric equation, that, in its general form, 
presents three regular singular points. According with the location of its 
singularities, the equation is known with different names (generalized 
Legendre equation. Gauss equation, ...). Confluence of two or all three 
of its singularities transforms the hypergeometric equation in others, 
with particular names (Kummer, Bessel, Weber, ...), encountered in 
the solution of a lot of problems. 

The next, in order of complexity, FDE is the Heun equation, that 
presents four regular singularities. Different kinds of confluence of these 
singularities produce the confluent, double confluent, biconfluent and 
triconfluent Heun equations. The book edited by Ronveaux (1995) con- 

© 2008 Kluwer Academic Publishers. Printed in the Netherlands. 



DCHeun_rv.tex; 7/11/2008; 16:29; p.l 



2 Abad e.a. 

Table I. Singularities of the different types of Heun equations 

equation rank rank 1 rank 2 rank 3 

Heun 0, 1, a (arbitrary ^ 0, 1, oo), oo 

confluent H. 0, 1 oo 

double confluent H. 0, oo 

biconfluent H. oo 

triconfluent H. oo 



tains a thorough revision of the results and open problems concerning 
those equations. More recently, Slavyanov and Lay (2000) have con- 
sidered again the class of Heun equations, in their study of Special 
Functions as solutions of differential equations with singular points. 
Confluence reduces the number of singularities, but increases the rank 
of some ones. Trivial changes of variable allow to locate the singularities 
at chosen points. We report in Table I the type of singularities and their 
positions when those equations are written in their usual conventional 
form . 

In almost all nontrivial applications, in which special functions arise, 
the problem at hand reduces to a boundary value one on an interval 
of the real axis whose interior points are ordinary ones, whereas one or 
both of its ends may be singular points of the differential equation. 
Usually, all but one of the parameters appearing in the coefficient 
functions of the equation are fixed. Then, the boundary conditions 
restrict (quantize) the possible values of the free parameter. In the 
neighbourhood of each end, a pair of independent solutions with a well 
known behaviour can be chosen. Such pair can serve as a basis to write, 
by linear combination, any other solution of the differential equation, in 
particular, that satisfying the boundary condition at the corresponding 
end. In order to examine the behaviour of this solution at the other end, 
one has to be able to write each one of the two solutions forming the 
basis at one end as linear combination of the basic solutions at the 
other end. In other words, one has to solve the so called "connection 
problem" for the ends. Articles by Kohno (1974), Naundorf (1976), 
Schafke and Schmidt (1980), and Schafke (1980, 1984), among others, 
deal with this problem. 

In a recent paper (Gomez and Sesma, 2007) we have suggested a 
different procedure to solve the connection problem for the cases of 
one of the ends being an ordinary or a regular singular point, and 
the other an irregular one. Examples of such cases are found in the 
triconfluent and biconfluent Heun equations, respectively. Our method 



DCHeun_rv.tex; 7/11/2008; 16:29; p. 2 



Double confluent Heun equation 



3 



is mostly inspired by the work of Naundorf, but presents some consid- 
erable advantages, from the numerical point of view. Our use of the 
Wronskians of the basic solutions to calculate the connection factors 
allows us to avoid having to sum slowly convergent double series like 
those encountered in the Naundorf's procedure. The purpose of this 
paper is to extend our method to the solution of the connection problem 
for two irregular singular points, the DCHE being a typical example of 
this case. 

We start by presenting, in the next Section, the DCHE and its basic 

solutions in the vicinity of its singular points. Our procedure to solve 
the connection problem and the corresponding algorithms are reported 
in Sections 3 and 4, respectively. An example, in Section 5, illustrates 
the method. Section 6 explains how to test the results of the only 
nontrivial algorithm. Finally, Section 7 contains a comment about the 
generality of the procedure. 



2. The double confluent Heun Equation 

A complete discussion of the different forms adopted by the DCHE, 
under appropriate transformations, has been done by Schmidt and Wolf 
(1995). We concentrate on what they call normal form (Eq. 1.1.10 of 
their article) 

, 2 

D^y + Biz)y = 0, D = z—, B{z)=Y,BpzP, (1) 

p=-2 

that, with the change of dependent variable 

wiz) = z'/^y{z) (2) 

and the notation 

Ao = Bo + l/A, Ap = Bp, pj^O, (3) 
can be written in the form 

p=-2 

free of first order derivatives and, therefore, preferable for our discussion 
because of the fact that the Wronskians of their solutions are indepen- 
dent of z. In what follows, we will assume that we are dealing with a 
non- degenerate case of the DCHE, that is, that 

A2 A_2 + 0. 
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The origin and the infinity axe the only singularities, both of rank 1, of 
the differential equation. We are interested in considering three pairs 
of independent solutions, namely: 

— Two multiplicative solutions (Arscott, 1987), wi and W2, that, ex- 
cept for particular sets of values of the parameters Ap in Eq. (4), 
have the form 

oo oo 

Wj = z''^ Cnjz'', being ^ |c„jp < oo, j = 1,2. (5) 

n=— oo n=— oo 

The indices Uj are not completely defined. They admit addition 
of any integer (with an adequate relabeling of the coefficients). To 
avoid ambiguities, we assume that 

< 1/2. 

— Two formal solutions, and W4, that have the nature of asymp- 
totic expansions for 2: — >■ 00, 

00 

Wk{z) =ex.p{akz)z^"' Y am,kz~'^, ao,k¥'0, k = 3,4. (6) 

m=0 

It is usual to say that these two expansions are associated to each 
other. 

— Two formal solutions, W5 and wq, asymptotic expansions for z — > 0, 
of the form 

00 

wiiz)=exp[Piz-')zP^ J2^rn,iz'", bo,l ^ 0, 1 = 5,6. (7) 

m=0 

Also these expansions are associated. 

The dctcrminaton of the index i' and the coefficients c„ of the multi- 
plicative solutions is rather laborious. By substitution of (5) in (4) one 
obtains the infinite set of homogeneous equations for the coefficients 

2 

{n+u){n-l+u)cn+ Y ApCn-p = 0, n =..., -1,0,1, ... , (8) 

p=-2 

that can be interpreted as a nonlinear eigenvalue problem, where the 
eigenvalue v must be such that 

oo 
n=— oo 
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In Section 4 we recall the Newton iterative method to solve that prob- 
lem. In general, two indices, I'l and 1^2, and two corresponding sets 
of coefficients, {c„^i} and {c„^2} are obtained, but for certain sets of 
values of the parameters Ap only one multiplicative solution appears. 
Any other independent solution must include powers of the variable 
multiplied by its logarithm. Such logarithmic solutions cannot corre- 
spond, usually, to the practical system that one tries to describe and 
are, therefore, to be discarded. We will assume, from now on, that 
the parameters Ap are such that the DCHE admits two independent 
multiplicative solutions. In what concerns the formal solutions, the ex- 
ponents a, 11, (5, and p and the coefficients and hm of the expansions 
(6) and (7) must be such that these expansions satisfy the differential 
equation. One obtains in this way 

"fe = V-M, Hk = -Ai/2ak, (10) 
Pi = ^^AT^, pi = l + A_,/2Pi, (11) 

for the exponents, and the recurrence relations 

2 

2akmam,k = {rn-iJ,k){rn-l-iJ,k) am-i,k + ^ A_pam-i-p,k, (12) 

p=0 

2 

2l3imbm,i = im-l+pi){m-2+pi)bm-i,i + ^Apbm-i-p,i, (13) 

p=0 

for the coefficients. The two independent solutions labeled by A; = 3, 4 
and those by Z = 5, 6 correspond to the two opposite values of the 
square root function in (10) and (11). 

3. The connection factors 

Any solution w of the DCHE can be written as a linear combination of 
the two multiplicative solutions, 

W = (iWi+(2W2. (14) 

Its behavior in the neighbourhood of the singular points can be imme- 
diately written if, besides the coefficients (i and (2, one knows the 
behaviour of the multiplicative solutions, that is, if one knows the 
connection factors T of their asymptotic expansions, 

Wj ~ Tj^3W3 + Tj^4^W4^, for z — ^ 00, j = l,2, (15) 
Wj ~ Tj^5W5 + Tj^QW6, for z^O, j = l,2. (16) 
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These connection factors are, obviously, numerical constants, but their 
values depend on the sector of the complex plane where z lies. This 
fact, known as "Stokes phenomenon" (Dingle, 1973), introduces a slight 
complication in the procedure. As it is well known, the connection 
factor multiplying each one of the asymptotic expansions in the right 
hand sides of (15) and (16) takes different values in the sectors of the 
complex z-plane separated by a Stokes ray of the associated expansion. 
On the ray, the value of the connection factor is the average of those 
two different ones. In the DCHE, both singular points have rank 1 and, 
consequently, for each one of the expansions (6) and (7), the sectors 
delimited by two contiguous Stokes rays of the same expansion have 
amplitude 27r. In fact, the Stokes rays for those expansions are given 

by 

arg z = — arg it 2n7r, n = 0, 1, 2, . . . , for lUjt, A; = 3, 4, (17) 



arg z = arg /3i it 2mr, n = 0, 1, 2, . . . , tor wi, Z = 5, 6 (18) 

To be specific, let us assume that the labels {3, 4} and {5, 6} are 
assigned in such a manner that 

either — tt < argQ;3 < — 7r/2 or 7r/2 < arg as < tt, (19) 

— 7r/2<argQ;4 < 7r/2, (20) 

either — tt < arg/^s < — 7r/2 or 7r/2 < arg (3^ < tt, (21) 

-7r/2 < arg/36 < 7r/2. (22) 

Then in the principal Riemann sheet, — vr < argz < tt, Tj^^ and Ij^4 
change their values as z crosses respectively the rays argz = — arg 04 
and arg 2; = — arg 0:3, whereas Tj^s and Tj^q change at arg 2; = arg/Jg 
and arg z = arg /Js , respectively. 



Now, we are going to present our procedure to calculate the con- 
nection factors. Let {r,s} denote either the pair {3,4} or {5,6}. Prom 
Eqs. (15) or (16), it is immediate to obtain 

_ W[lVj,Ws] _ W[Wj,Wr] , . 

where W[f,g] represents the Wronskian of the functions / and g, 

W[f,g]{z) ^ f{z)g'{z) - f'{z)g{z). (24) 

There is no difficulty in using the asymptotic expansions (6) and (7) to 
calculate the denominators in (23). In fact, 

W[w3,W4\ = ->V[w4,ti;3] = 20400,300,4, (25) 
yV[w5,W(i] = - W[wg, W5] = 2/35 60,5 bo,6. (26) 
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The numerators, instead, require a more sophisticated method. Direct 
evaluation would give asymptotic expansions, containing both positive 
and negative powers of the variable, which result to be inadequate to 
obtain precise values of those numerators. But we can benefit from the 
fact that, whenever one considers not merely a ray, but a sector of the 
complex plane, an asymptotic power series defines unambiguously an 
analytic function. Our idea is to write, for each one of the connection 
factors Tj^t {t = r,s), two asymptotic expansions of the same function, 
one of the expansions containing, as a common multiplicative constant, 
the Wronskian in the numerator of the expression of Tj^f Comparison 
of coefficients of equal powers of the variable in the two expansions 
gives immediately that numerator. With that purpose, we introduce 
some auxiliary functions and asymptotic expansions, namely 

Uj{z) = e-"*^/2 ^.^ = j = 1, 2, t = 3, 4, (27) 

Vj{z) = e-^*/2^ Wj, vt{z) = e-^*/2^ j = 1, 2, t = 5, 6, (28) 
whose Wronskians obey the relations 

W[uj,ut] = e-"'' W[wj,wt], J = 1,2, t = 3, 4, (29) 
W[vj,vt] = e-^'-'~'w[w,,wt], j = l,2, t = 5,6. (30) 

Direct calculation of the left hand sides of these equations gives the 
expansions 

oo 

W[uj,ut] = E 7n,,-t^"^"^+''S j = l,2, t = 3,4, (31) 

n=— oo 
oo 

mvj,vt] = E ^"J.t^"^"'^"*' ■?' = 1'2, * = 5,6, (32) 

n=— oo 

with coefficients 

oo 

1n,j,t = E (^m,t[oit Cn+m,j - in+2m+l + l'j - fJ,t)Cn+m+l,j], (33) 

m=0 

oo 

^n,j,t = E bm,t[-PtCn-m+2,j - (n-2m+ 1 + -/9t)Cn-m+ ,(34) 

m=0 

Now, in order to compare coefficients of equal powers in the two sides, 
we need to write expansions of the right hand sides of (29) and (30) in 
powers of z with the same exponents as in the left hand sides. This can 
be done by using the Heaviside's exponential series (Hardy, 1949) 

oo fn+5 

^ J2 Y(n + l + d) ' |argC|<7r, 5 arbitrary, (35) 

n=— oo \ ' ' f 
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with ^ replaced respectively by —atz and —I3tz~\ and taking 5 equal 
to lyj+Hk a-iid —(uj + pt). This kind of representation of the exponential 
function has already been used by Naundorf in his solution of the 
connection problem. Notice that the restriction | arg^| < vr prevents the 
use of such expansions if argz = — argaj or argz = aigPf But these 
values of arg z correspond precisely to the rays at which the connection 
factor to be calculated changes. As we have already said, the value 
assigned to the connection factor in that situation is the average of 
its values in the sectors separated by that ray, where the restriction 
I arg^l < TT is fulfilled. Following that procedure, one obtains, for z out 
of the Stokes ray of wt in the principal Riemann sheet, 

= ^t-^lt'iZ'^ ^nj,t, J = 1, 2, t = 3, 4, (36) 

W[wj,wt] = ^|^^^)~'^.~?^ -nj,t, J = 1,2, t = 5,6, (37) 

where the minus sign in front of at and Pt is to be interpreted as e*'^ 
or e~"^ so as to have | arg(— at2;)| < vr and | arg(— /3iZ~-^)| < tt. If z lies 
on the Stokes ray of Wt, one has, instead of (36) and (37), 

W[w„wt] = (-l)-cos[7r(^.,+/.,)]^^^^^±i±^±^7n,,,i, (38) 

at 

W[wj,wt] = (-1)" cos[7r(^,- +p,)] ^^""^nZX^'^ (39) 

Pt 



4. The algorithms 

Wc present in this Section algorithms that facilitate implementation 
of the procedure sketched above. Its successive steps are considered 
separately in different subsections 

4.1. Multiplicative solutions 

As it has been already said, we face in Eqs. (8) and (9) a nonlinear 
eigenvalue problem. Algorithms to solve finite order problems of this 
kind have been discussed by Ruhe (1973). The condition (9) implies 

liin \cn\ = 0, (40) 

n— >±oo 

that makes possible to reduce, by truncation, our problem (8) to one 
with n going from —M to N, both M and A'' being positive and 
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sufficiently large to ensure that the solution of the truncated problem 
approximates that of the original one. Each step of the Newton iteration 
method consists in moving from an approximate solution, c4^^}, 
to another one, {v'^^~^^\ ch~^^^ } , by solving the system of equations 

(2n-l+2z.«)cW(^(^+i) - ^.W) + (n+z.«)(n-l+i.«)c(:+i) 

2 

+ J2 Ap^n-p =0, n = -M, . . . , -1,0, 1, . . . ,iV, (41) 

p=-2 

N 

E <^*4^'^ = 1' (42) 

n=-M 

that results, by linearization, from (8) and from the truncated normal- 
ization condition 

N 

n=-M 
(i) 

Needless to say, the values of c^n with m < —M oi m > N entering in 
some of the equations (41) should be taken equal to zero, in accordance 
with the truncation done. As usual, the iteration process is stopped 
when the difference between consecutive solutions is satisfactory. Then 
the process is repeated with larger and larger values of M and N, to 
obtain a stable solution within the required precision. 

The iteration process just described needs initial values {v^^\cn^} 
not far from the true solution. The two different values of v can be 
obtained from the two eigenvalues 

\j = exp(2i7ri/j) (43) 

of the circuit matrix C (Wasow, 2002) for the singular point at z = 0. 
The entries of that matrix can be computed by numerically integrating 
the DCHE on the unit circle, from z = exp(O) to z = exp(2i7r), for two 
independent sets of initial values. If we consider two solutions, Wa{z) 
and wi,{z), obeying, for instance, the conditions 

waie"^) = 1, w',ie'>) = 0, 
^6(eO) = 0, w;Ue°) = 1, 



then 
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and 

1 



In 
2i7r 



7: Cii + C22 ± \] {C^l-C22f + ^Ci2C2l 



(44) 



The two signs in front of the square root produce two different values for 
z/, unless the parameters Ap in the DCHE be such that (Cn — C22)^ + 
4C12C21 = 0, in which case only one multiplicative solution appears, 
any other independent solution containing logarithmic terms. The am- 
biguity in the real part of v due to the multivaluedness of the logarithm 
in the right hand side of (44) is eliminated by the restriction \^v\ < 1/2 
assumed above. Notice that 

Ai A2 = det C = W[wa, Wb] = 1 

and, therefore, 

I'l + 1^2 = (mod 1). 

This may serve as a test for the integration of the DCHE on the unit 
circle. 

Although Eq. (44) is exact, the Cmn are obtained numerically and 
the resulting values of u may only be considered as starting values, Vj^^ 

{j = 1, 2), for the Newton iteration process. As starting coefficients 
one may use the solutions of the homogeneous system 

(n+.f)(n-l+.f)ciVE^.^i\.=0 

p=-2 

^ = _M,... ,-1,0,1,..., AT, j = l,2, (45) 
with the already mentioned truncated normalization condition 

N 



M |2 

n=—M 



4.2. Formal solutions 



The exponents of the formal expansions (6) and (7) are given in (10) 
and (11). The coefficients must obey the recurrence relations (12) and 
(13) that are but third order difference equations. The Perron-Kreuser 
theorem (Perron, 1959) predicts for each one of them a unique (save for 
multiplication by a constant) dominant solution that can be obtained 
starting, for instance, with 

ao,k = l, bo,i = l, A; = 3,4, / = 5,6, (47) 
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and using the recurrence relations as they appear in (12) and (13). 
To avoid overflows, it may be convenient to deal with the quotients of 
successive coefficients, that satisfy 

"^^km-^^ = {m-Hk){m-l-Hk) + Ao 

"m.-l,/c Qm-l,fe Qm-2,fe ' \^°) 

°'m-2,k "m-2,fc "m-3,fc 



2Am5^ = {m-l+pi){m-2+pi)+AQ 



Om-l,l Om-l,l Om-2,l ^ ' 

bm-2,l bm-2,l bm-3,i 



4.3. Connection factors 

For the computation of the connection factors one should make use of 
Eqs. (23), (25), (26), (36) and (37), complemented by (33) and (34) 
conveniently truncated. The integer n in Eqs. (33), (34), (36) and (37) 
may be chosen at will, with the only restriction of being positive and 
satisfying 

2 

{n + iy){n + i^-l)> ^ \Ap\. (50) 

p=-2 

Use of different values of n may serve as a test of the procedure. The 

sums in (33) and (34), truncated and written in terms of quotients of 
successive coefficients, am,k/ o.m-\,k a^id &m,//^m-i,/) niay be computed 
in nested form. 



5. An example 

A particular case of DCHE on the positive real semiaxis, z G [0, +oo), 
is the Schrodinger equation 

-|^(^-^i^M)+v-Mij(,., = £iiH. (51) 

for the reduced radial wave function R{r) of a particle of mass m, 
angular momentum Ih and energy E = A2ft^ /2mrQ in a spherically 
symmetric potential 

Vir) = -f(^ + ^ + ^!^±^ + ^]. (52) 
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In fact, by using the variable z and the wave function w given by 

z = t/tq and w{z) = R{r), (53) 

the Schrodingcr equation (51) adopts the form (4). The potential (52), 
with some restrictions on the values of the parameters, belongs to a 
class of quasi-exactly-solvable ones (Turbiner, 1988). In fact, for the 
particular set of parameters 

A_2 = - 1, A_i = 4/5, Ao = 31/25, = 3/5, (54) 

it presents an / = bound state of energy E = — (1/4) /i^/2mro 
(Ozgelik and §im§ck, 1991). In other words, Eq. (4), with the values of 
the Ap given by (54) and by A2 = —1/4, has a normalizable solution 
on the positive real semiaxis. To illustrate the procedure described in 
the preceding Sections, we have applied it to find global solutions, for 
z G [0, -l-oo) , of (4) with fixed values (54) of the potential parameters 
and several values 

A2 = - 1/10, - 1/5, - 1/4, - 3/10, - 2/5, (55) 

of the energy parameter. The results are shown in Tables II to V. We 
report the output of our double precision FORTRAN codes, but, due 
to roundoff errors, we do not claim that all the digits reproduced arc 
correct. In fact, entries that in the tables appear as having modulus less 
than lO"-*^^ should be exactly equal to zero. Table II shows the index 
vi of the multiplicative sohition wi. Of course, the index V2 of W2 is 
the opposite, 1^2 = —i^i- The connection factors of both multiplicative 
solutions with the formal ones are listed in Table III. These factors 
depend on the normalization adopted for the different solutions. We 
have taken 

coj = 1, i = l,2, 

for the multiplicative solutions, whereas the normalization of the formal 
solutions is determined by 

ao,k = 1, = 3, 4, 6o,i = 1, ^ = 5, 6. 

Table IV gives the coefficients Ci and (2 of a linear combination of the 
multiplicative solutions, 

wreg = Ciw'i + C2W^2, 

well behaved near the origin on the positive real semiaxis, that is, such 
that 

wxeg{z) ^ W5{z) as z ^ 0"*". 
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Table II. Indices of the multiplicative solutions of the DCHE for the sets of param- 
eters given in (54) and (55). The two indices have opposite values; so the real and 
imaginary parts of only one of them are shown. 



A2 






-1/10 


-.500000000000E+00 


-.703150555392E+00 


-1/5 


-.500000000000E+00 


-.531738153810E+00 


-1/4 


-.400000000000E+00 


.102867462041E-26 


-3/10 


-.262974969075E-15 


.509507933497E+00 


-2/5 


.181198729462E-16 


.688682990633E+00 



Finally, Table V shows the connection factors giving the behaviour of 
■ii^reg at infinity on the positive real semiaxis, 

wveg{z) ~ Treg,3 Wi{z) + Treg,4 Wi{z) as z ^ +oo. 



We have already mentioned that the DCHE (4) with parameters as 
given by (54) and by ^2 = — 1/4 possesses a solution, 



w{z) = z^l^ exp(-z-^ - z/2), 



(56) 



normalizable on the positive real semiaxis. This fact allows to know the 
exact values of the connection factors Ti^a and Ti^s of the multiplicative 
solution 



oo 
n=— oo 



n-2/5 



Writing (56) in the form 



(57) 



n+3/5 



(58) 



with 



Cn = < 



(-If £^7 ^ forn<0, 

m'.im — n)\ 

(-irE-77— ^ forn>0. 



m\{m + n)! 
comparison with (57) gives immediately 

T^ifi = Ti,5 = co,i/c_i 
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Table III. Connection factors of the multiplicative solutions, wi and w-z, with the formal 
solutions at infinity, W3 and W4, and at the origin, W5 and Wg, for the DCHE (4) with 
parameters (54) and (55). 


A2 






5Rri,4 


3Ti.4 


-0.1 


-.130166356702E+01 


.218375080422E+01 


537554017411E+00 


.314875723054E+00 


-0.2 


-.204239771333E+01 


.362437608911E+00 


462340626214B-01 


.190707221682E+00 


-0.25 


- .786334477859E+00 


. 770822392694B - 13 


227716282126B-13 


.938023374790B-28 


— 0.3 






iUozi8 i^8;>9ic. + UU 




-0.4 


.108556195947E + 01 


.199397670471E+01 


125046997516E+00 


.673923636677E-01 


A2 


KT2,3 




5RT2,4 


ar2,4 


-0.1 


-.179473952568E+00 


.141771950373E+01 - 


.347068293144E+00 


-.466342783606E-01 


-0.2 


.104106884981E+01 


.818532573406E+00 


834077925914E-01 


.934810009512E-01 


-0.25 


.491548162537E+00 


.488371898613E-13 


289740699463B-01 


-.293358750439B-14 


— 0.3 


1 rt'7C /I c c n*7C c 1 nn 
. lU704D0ooy tOtj-j-UU 




iUDZi8 / 0809iil + UU 


— . iU8 /13b8z:J4Jllj — Ui 


-0.4 


.108556195947E+01 


-.199397670471E+01 


125046997516E+00 


-.673923636677E-01 


A2 


»Tl,5 




KTi,6 


3^1,6 


-0.1 


-.3S0619618820E+01 


-.780802661791E+00 


283770014388E-01 


- . 136520514662E+00 


-0.2 


-.218278236288E+01 


-.183368690581E+00 


349930523633E-02 


- .829409172034B-01 


-0.25 


- . 786334477859E+00 


- . 770774921662E- 13 


581790865708B-15 


.371252705024B-27 


-0.3 


- . 185903592530E+00 


.135138118813E+01 


590676874067B-01 


.526872360606B-02 


-0.4 


-.379273027228B-01 


.219921373757E+01 


927424524732E-01 


.159396234117B-03 


A2 










-0.1 


.182308938715E+01 


.8679382D8517E+00 


346916923655E-01 


-.702844924334E-01 


-0.2 


.128432288745E+01 


.553414307656E+00 


229815763899E-01 


-.477580284950E-01 


-0.25 


.493413857062E+00 


-.477076860664B-13 


144870349731B-01 


- . 142012480382B- 14 


-0.3 


- . 185903592530E+00 


-.135138118813E+01 


590575874067B-01 


-.626872360505B-02 


-0.4 


- .379273027228E-01 


-.219921373757E+01 


927424524732B-01 


- . 169396234117B-03 



and, having chosen co,i = 1, 

Ti,3 = Ti,5 = — 

The values of 3fiTi,3 and sRTi^g 
ported in Table III coincide w 
digits. 



oo n-m ^ 

y — = — ■ 

J^^m\{m + l)\_ 

obtained with our procedure and re- 
th the exact values up to all shown 
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Table IV. Coefficients of the linear combination of multiplicative solutions resulting in a 
solution wreg = Ci^i + C2W2 normalizable, on the positive real semiaxis, at « = 0. The 
parameters of the DCHE are the same as in Tables II and III. 



A2 






KC2 


aC2 


-0.1 


-.136295750328E+00 


.283302821758E-01 


.222074702056E+00 


-.109613756594E+00 


-0.2 


-.2282B6434722E+00 


.963021687162E-02 


.322449462131E+00 


- . 15516B46B3B8E+00 


-0.25 


-.12717234B631E+01 


.124665928872E-12 


.610716706831E-13 


.325901912099E-26 


-0.3 


-.325468718708E-01 


-.365514500226E+00 


-.325468718708E-01 


.365514500226E+00 


-0.4 


-.390741170002E-03 


-.227347243108E+00 


-.390741170002E-03 


.227347243108E+00 



Tabic V. Connection factors of the "regular" solution at the origin with the formal 
solutions W3 and W4 for a ^ 00 along the positive real semiax;is. The parameters in 
the DCHE are the same as in the preceding Tables. 



A2 




^ ?reg,3 


3^ -^reg,4 






-0.1 


.231089872113E+00 


-.486000928307B-12 


- . 164373692468B+00 


.317038062470E-13 


-0.2 


.926400136830E+00 


.898669004980E-12 


-.24779B48019BE- 


■01 


.B60940183192E-13 


-0.25 


.1()()()()()000000E+01 


-.196048546485E-12 


-.304389678921E- 


13 


.162,'523413961E-26 


-0.3 


.100006470514E+01 


.OOOOOOOOOOOOE+00 


.148614535378E- 


01 


.OOOOOOOOOOOOE+00 


-0.4 


.906801865774E+00 


.OOOOOOOOOOOOE+00 


.307406681930E- 


01 


.OOOOOOOOOOOOE+00 



6. A test of the multiplicative solutions 

Wc have described, in subsection 4.1, an algorithm to obtain the mul- 
tiphcative solutions, wi and -102, mentioned in Eq. (5). A test of the 
correct implementation of the algorithm could be to try to reproduce 
the particular solution of the DCHE provided by the Maple 11 system. 
The statement 

HeunD(Q;, /3, 7, 5, t) 

gives the value at t {\t\ < 1) of a function y{t) that obeys the DCHE, 
written in the Jaffe-Lay form (Slavyanov and Lay, 2000), 

<fy{t) _ a + 2t + at^ - 2t^ dy 6 + {2a + j)t + (5t^ _ 

dt^ (t2 - 1)2 dt + (t2 - 1)3 y- ^^^> 

with singular points located at i = — 1 and i = 1, and satisfies the 
boundary conditions 

2/(0) = 1, y'(0) = 0. (60) 
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(We adopt the notation of the Maple manual. Needless to say, the 
parameters represented by the symbols a and /? along this Section 
are not related to the exponents ct^, k = 3,4, and / = 5, 6 of 
the formal solutions (6) and (7).) The changes of independent and 
dependent variables 

t = ^ and y{t) = z'^/^ g - ^) ) w{z) (61) 

transform Eq. (59) in Eq. (4) with parameters 

«2 J-/3-6 8-«^ + 4(/3-(5) 
^-^ = "64' = 16 ' = 32 ' 

16 ' 64 ^ 

Obviously, w{z) in (61) can be written as a linear combination 

w{z) = ^1 ^1(2;) + ^2 W2{z) (63) 

of the two multiplicative solutions, the coefficients ^1 and ^2 depending 
on the parameters a, (5, 7 and 5. Then, the procedure to check the 
algorithm used to obtain the multiplicative solutions would run along 
the following steps: 

— Choose a set of values for the parameters a, /3, 7 and 5. 

— Compute the corresponding values of Ap, p = — 2, — 1, . . . , 2, by 
means of the relations (62). 

— Use the algorithm described in Subsection 4.1 to obtain the two 
multiplicative solutions. 

— Determine the coefficients ^1 and ^2 in (63) by requiring that the 
second of Eqs. (61) and that resulting by deriving it with respect 
to t be satisfied for t = Q and z = 1. In other words, solve the 
system 

wi{z = l)^i+W2{z = l)i2 = 1, (64) 
w'i{z = l)ii+ w'^iz = 1) 6 = - 1/2. (65) 

— Choose an arbitrary value of t inside the unit circle and compute 
the corresponding value of z, according to the first of Eqs. (61). 

— Check the fulfilment of the second of Eqs. (61) for those val- 
ues of t and z, the left hand side being computed by the above 
mentioned Maple statement and the right hand side by using the 
multiplicative solutions and the coefficients ^1 and ^2- 
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Table VI. Coefficients Cn,j of the multiplicative solutions, Wj{z) = z"^ X^^_oo '^^i 
j = 1,2, normalized in such a way that co,i = Co, 2 = 1- Only the most relevant 
coefficients axe shown. The indices are vi = — V2 = 0.346120772343. 



m 


C— m,l 




Cm,l 




C-m,2 




Cm, 2 




1 


. 167014439267E+00 


.190993398030E+01 


- . 107279310787E+00 


.382083714867E+00 


2 


.481376140064B- 


01 


-.601406710996B- 


02 


.396618603946B- 


01 


.343615502071E+00 


3 


.266845662796B- 


02 


.651050367942B- 


01 


-.297089364186B- 


02 


.126405802915B- 


01 


4 


.670742341737E- 


03 


-.688158578777B- 


03 


.48411ir)13331E- 


03 


.916616846207E- 


02 


5 


.151109427662E- 


04 


.719872336102B- 


03 


-.276617916112E- 


04 


.120052304745E- 


03 


6 


.436467584453B- 


05 


-.788297814632B- 


05 


.284956022379B- 


05 


.881920361011B- 


04 


7 


.424169499646B- 


07 


.392398781019B- 


05 


-.131313020620B- 


06 


.507159044452B- 


06 


8 


.163574045519E- 


07 


-.405522033871E- 


07 


.982771078139E- 


08 


.435958940659E- 


06 


9 


.537668556535E- 


10 


.127234781325E- 


07 


-.379605967044E- 


09 


.109230687414E- 


08 


10 


.397572715764B- 


10 


-.122029907678B- 


09 


.222381579709B- 


10 


.13ia80043658B- 


08 


11 


-.119341203300B- 


13 


.273442688328B- 


10 


-.739683873667B- 


12 


.10S746798402B- 


11 


12 


.676604266157B- 


13 


-.243166143611B- 


12 


.366327272482B- 


13 


.264897017672B- 


11 


13 


-.164661027834B- 


15 


.417999894189B- 


13 


-.103816479336B- 


14 


-.464a27978928B- 


15 


14 


.860622285001B- 


16 


-.346793838238B- 


15 


.422141078690B- 


16 


.384446086066B- 


14 


15 


-.319253539778E- 


18 


.477715538676E- 


16 


-.109995830505E- 


17 


-.29969646()688E- 


17 


16 


.821887767466E- 


19 


-.369227215875E- 


18 


.387438635716E- 


19 


.420067427815E- 


17 


17 


-.376967703210B- 


21 


.423682849020B- 


19 


-.911209342016B- 


21 


-.490269012114B- 


20 


18 


.629190033709B- 


22 


-.307216713853B- 


21 


.282929069113B- 


22 


.358011576237B- 


20 


19 


-.321864802770E- 


24 


.29986644121 IE - 


22 


-.606469439157E- 


24 


-.509135489559E- 


23 


20 


.390955943808E- 


25 


-.204922402639E- 


24 


.168292193160E- 


25 


.244636867112E- 


23 



Let us choose, for instance, the parameters 

a = 4, /? = -3, 7 = 2, (5 = -l, 

for which the corresponding parameters in Eq. (4) are 

^_2 = -l/4, ^-1 = 3/8, Ao = -l/2, ^1 = 1/8, ^2 = -1/4, 

The multiplicative solutions, obtained with a double precision FOR- 
TRAN code, are then those given in Table VI, and the coefficients in 
(63) giving the w{z) corresponding to the Maple particular solution 
turn out to be 

^1 = - 0.48013092979925, ,^2 = 1.5087428376316. 
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7. Final comment 



The procedure presented in the preceding Sections referred to global 
solutions in the interval between two adjacent irregular singular points, 
of rank 1, of a second order differential equation like the DCHE. It 
can be trivially extended, however, to the case of two irregular singular 
points of arbitrary rank. The only changes due to the larger rank of the 
singularities are in the evaluation of the Wronskians of the pairs formed 
by each one of the multiplicative solutions and each one of the formal 
solutions. The steps are analogous to those detailed in a paper (Gomez 
and Sesma, 2007) dealing with the connection problem for an ordinary 
or regular singular point and an irregular one of arbitrary rank. Obvi- 
ously, the series solutions around the ordinary or regular singular point 
in that reference should be replaced by the multiplicative solutions of 
the problem at hand. Then, having computed the connection factors of 
each multiplicative solution with the formal solutions at each singular 
point, the procedure is the same as described above. 
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